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A measurement of the Lamb shift in muonic hydrogen yields a charge radius of the proton that 
is smaller than the CODATA value by about 5 standard deviations. We explore the possibility that 
new scalar, pseudoscalar, vector, and tensor flavor-conserving nonuniversal interactions may be re- 
sponsible for the discrepancy. We consider exotic particles that among leptons, couple preferentially 
to muons, and mediate an attractive nucleon-muon interaction. We find that the many constraints 
from low energy data disfavor new spin-0, spin-1 and spin-2 particles as an explanation. 



Lamb shift. The success of quantum electrodynam- 
ics (QED) is apparent in the explanation of the Lamb 
shift |l[ which is the observation that the 2Si/2 state of 
hydrogen is higher than the 2P 1 / 2 state by about 1 GHzQ 
Precision measurements in atomic spectra have tested 
bound-state QED to the extent that the charge distribu- 
tion of the proton needs to be taken into account Q • The 
root-mean-square charge radius of the proton compiled 
by CODATA from the spectroscopy of atomic hydrogen 
and electron-proton scattering is [3| 
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0.8768 ±0.0069 fin, 



(1) 



provided there are no new long-range e-p interactions [J] . 
It has been a long-held goal to measure the correspond- 
ing Lamb shift in muonic hydrogen which is even more 
sensitive to the structure of the proton due to its smaller 
Bohr radius (am,,)" 1 , (where a ~ 1/137 is the electro- 
magnetic fine structure constant and m M ~ 105 MeV). 
Recently, the 2P^^ 2 — > 2S[^ Lamb shift in muonic hy- 
drogen was measured to be [5j , 

AE = E(2P£ff)-E{2Slj?) = 206.2949±0.0032meV, 

(2) 



while the predicted value is [H, 
AE = 209.9779(49) - 5.2262 (r 2 ) + 0.0347 (r 2 ) 
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(3) 



1 The dominant contributions to the Lamb shift arise from vac- 
uum polarization and the vertex charge form factor of the lep- 
ton. Vacuum polarization contributes negatively to AE = 
E(2S\/2) — E(2Pi/2) since more of the lepton's bare charge is 
revealed for the S state (than the P state) due to its greater 
overlap with the nucleus. On the other hand, the vertex charge 
form factor is related to the zitterbewegung of the lepton which 
causes the effective Coulomb potential to be smeared out and less 
attractive. The effect is greater for the S state, so that the con- 
tribution to AE is positive. The latter contribution is dominant 
for ordinary hydrogen but plays a minor role in muonic hydrogen 
because of the smaller Compton wavelength of the muon. Con- 
sequently, AE is positive in ordinary hydrogen and negative in 
muonic hydrogen. 



where radii are in fm and energy in meV, and the number 
in parenthesis indicates the la uncertainty of the last two 
decimal places of the given number. (Note that AE is 
defined to be positive.) Eqs. ([2]) and (|3|) yield the order 
of magnitude more precise result [5j , 



(r 2 ) 1/2 = 0.84184 ± 0.00067 f m , 



(4) 



which is smaller than the CODATA value by about 5a. 
A partial resolution of the discrepancy may be found in 
a correlation between (r 2 ) and the r 3 -dependent third 
Zemach moment (since they contribute to the Lamb shift 
with opposite signs) and perhaps unreliable extractions 
of these from electron- proton scattering data @ . Never- 
theless, a 4cr difference remains. The possibility that the 
4% difference is a hint of a new gauge interaction with a 
natural scale am^ has been entertained in Ref. 0. 

In this Letter, we postulate the existence of a new in- 
teraction between muons and nucleons, and study its na- 
ture, bearing in mind the many experimental constraints. 
The interaction must be attractive since AE measured in 
muonic hydrogen is larger than expected, signaling that 
the 2S\/2 state is subject to a stronger attraction than 
electromagnetic. 

Scalar and spin-2 boson exchanges produce an attrac- 
tive potential, giving positive contributions to AE. Pseu- 
doscalar boson exchange is a derivative interaction in- 
volving the spins and velocities of the lepton and the 
nucleus, which becomes insignificant in the nonrelativis- 
tic limit, and irrelevant to the Lamb shift. A boson with 
both scalar and pseudoscalar couplings violates CP con- 
servation. Such a scenario faces strong constraints from 
electric dipole moment measurements of leptons and nu- 
cleons, and is disfavored [T^. Vector boson exchange 
(like photon exchange) can produce an attractive poten- 
tial if the quantum numbers associated with the lepton 
and the nucleus are opposite in sign. Then Lamb shift 
phenomenology is like that of scalar exchange. Axial- 
vector exchange couples the spins of the lepton and the 
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nucleus in the nonrelativistic limit (with an effective po- 
tential —a x (cr Pi -(T p )e~ rnxr /r [H|), an d affects the hyper- 
fine structure (but not the Lamb shift) so that the cor- 
rection to the hyperfine splitting between the 2Py= 2 and 

2Sf^ 1 levels for to v 3> am u is a^^f- 1 + 1 °( m x/( Q;m ^ 

1/2 a p a 20 (Hm x /(dm r )) 4 ' 

in the notation defined below. However, since the axial- 
vector current is not conserved, the propagator gives a 
very singular contribution for m x < am^, which is un- 
physical. Absent a well-defined model, we do not consider 
the axial-vector case any further. 

Suppose the interaction between fermions / and \ is 
given by CV ' ffx, where S, V,T denote scalar, vector 
and tensor respectively, and / can be a muon ^ or 
a nucleon n; we assume isospin is conserved. Through- 
out, we take the couplings C to be real and positive. In 
the nonrelativistic limit, the muon-nucleon interaction is 
given by the Yukawa-type potential, 



AV(r) 



(5) 



where a x = C^' V ' T C^' V ' T / (4tt) and m x is the mass of the 
particle x- Physical systems in which tensor interactions 
(by which we mean spin-2 exchange) are governed by 



a Yukawa potential allow the identification, Cj 



Cf. 



This will be valid for all our constraints except those that 
involve the anomalous magnetic moment of the muon 
= (oa — 2)/2. The correction to the muonic Lamb 

shift is Tain 



S(AE) 
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where m r is the reduced mass of the muon-proton system, 
and its use is numerically important as m r is smaller 
than by more than 10%. The green shaded region in 
Fig.[T]shows the 95% C. L. region that accommodates the 
difference between Eqs. ([TJ and ((4]). We do not consider 
m x > 10 GeV since the required a x becomes larger than 
20a, in the nonperturbative regime. 

Upsilon decay. For scalar \ m the mass range 
2m^ — 9.3 GeV, the nonobservance of radiative decays of 
the T(25) and T(3S) resonances: T — » 7Y, x ~^ A f+ 
strongly constrains the T — \ coupling [13| , which we ex- 
pect to be no smaller than C„ ; for a Higgs-like x, the cou- 
pling is naturally C(mb/m„) x C„ , while for a universal 
interaction, it should be 0(C„). We conservatively take 
the T — \ coupling to be Cf . In obvious notation [l4[ , 



BF(T ->■ fi+fJ,-) 4?ra 



1-4 ■ (7) 



Under our assumption that the branching fraction of x ~ ► 
fi + fi~ is unity, the 90% C. L. upper limit on ranges 
from (0.94 — 9.4) ■ 10 -3 13J, where the lower end of the 
range corresponds to smaller m x . The values of a x /a 



53 





• // 

* // 


Vector 


* // 
» / 
* // 


Scalar 

Tensor 


' // * 

* It * 

* It * 

#/ 
/// 

// 
ft 




// 

V 

r 






"—" /// 





m x (MeV) 



FIG. 1. The 95% C. L. range of a x /a required to reproduce 
the muonic Lamb shift is indicated by the green shaded re- 
gion. The black solid, red dashed and blue dot-dashed lines 
are the upper limits for vector, scalar and spin-2 particles, 
respectively, from a combination of n— 208 Pb scattering data 
and the anomalous magnetic moment of the muon. The black 
dotted curve is the upper bound obtained from atomic X-ray 
transitions. All bounds are at the 95% C. L. 



needed to explain the muonic Lamb shift constrain to 
lie above 0(1), O(10) and 0(100) for m x - 2m M , 1 GeV 
and 9 GeV, respectively, couplings which are too large. 

A vector x can mediate leptonic decays of spin-1 
quarkonia. Since the only lepton that x couples to is 
the muon, one expects nonuniversality in leptonic decays. 
(6) For T(1S) decays, R T ^ = T TT /T^ = 1.005 ± 0.013 ± 
0.022 pH, whereas the SM expectation is 0.992 [lij. The 
inclusion of \ modifies the SM value of i? TM by a factor, 



1± 



uQt 



(m x /mx) 



1 — {m x /mfY 



(8) 



where the + (— ) sign corresponds to destructive (con- 
structive) interference and Qb is the electric charge of 
the b quark. For m x < 1 GeV, a conservative 95% C. L. 
(one-sided) upper bound on a x /a (assuming the SM and 
X contributions destructively interfere) is 8.8 x 10~ 3 . In 
the range 1 GeV < m x < mx , the upper bound becomes 
even more stringent, falling monotonically with m x . The 
mass of a vector x is restricted to be less than about 
230 MeV in order to explain the muonic Lamb shift. 

Neutron scattering. Very precise neutron scattering 
experiments on heavy nuclei in the keV regime have been 
performed to study the electric polarizability of the neu- 
tron. The goal is to measure interference effects between 
the nuclear potential and the r -4 potential produced by 
electric polarizability. One can then see that a Yukawa 
potential (Cf > V -T) 2 e ~ m * r / (Airr) may also be probed 
by such experiments; the minus and plus signs apply to 
scalar/tensor and vector interactions, respectively. Strin- 
gent bounds are obtainable because the p-wave amplitude 
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due to the short range strong interaction depends linearly 
on energy and differs markedly from that due to the new 
longer range interaction. A n— 208 Pb scattering experi- 
ment in the neutron energy range 1 to 26 keV mea- 
sured the differential cross section (under the assumption 
that the scattering amplitude can be expanded in s and 
p waves) to be 



da/dVt = <T (1 +uEcos6) /(47r) , 



(9) 



with s/ao/iTT^ lOfmandw = (1.91±0.42)-1CT 3 keV -1 . 
The measured values are in line with expectations so that 
the Yukawa potential contribution ought to be subdom- 
inant. Denoting the strong interaction contribution to 
to by u> s , and the contribution of the new interaction by 
Aw, clearly, uj = lo s + Acu with [l8| 



Aoj = =f 



16 (C^ T ) 2 Ami 



4tt 



(10) 



in the Born approximation (not valid for m x ^0.1 MeV), 
and m n is the neutron mass and A is the atomic mass 
number. For scalar/tensor exchange, it is possible that a 
cancellation between uj s and Au; produces the experimen- 
tal result. However, this cannot be the case for a vector 
X- A conservative 95% C. L. (one-sided) upper limit can 
be obtained by requiring that Aw < 2.6 x 10~ 3 keV, i.e., 

C v n < (m x /206) 2 , (11) 

with m x in MeV. It is the shaded region of Fig. [2] 

While reliable bounds for a scalar/tensor \ are n °t 
extractable from the differential cross section, the total 
cross section measured between 10 eV and 10 keV 3, 2(| 
may be employed with confidence. The energy depen- 



dence of the 



208 



Pb cross section for neutron energies 



below 10 keV can be parameterized by 

Cl(fc) = cr(O) + O- 2 fc 2 + 0(fc 4 ) 



(12) 



where k = 2.1968 x 10~ 4 VEA/(A+ 1) is the wave vector 
of the incoming neutron (with k in fm and E in eV). 
The cross section in the limit of vanishing momentum 
transfer ct(0) is directly related to the scattering length, 
and (72 gives the effective range of the potential. The 
0(k) contribution to a(k) arises from the electric field 
of the nuclear charge distribution and is negligible. The 
measured values cr(0) = 12.40 ± 0.02 b and a 2 = -448 ± 
3 b fm 2 give a 95% C. L. bound on Cf T [U that is 
almost identical to Eq. (fTTj) in the mass range of interest 
(and is not shown separately in Fig. [2]), but without the 
ambiguity from the cancellation mentioned above. 

Muon anomalous magnetic moment. We now 
consider the independent constraint on from a^. In 
fact, since the experimental value of is above the 
SM expectation by more than three standard deviations: 
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FIG. 2. The curves are 95% C. L. upper bounds on the muonic 
couplings C^, CY Cj from Aa M . The green shading marks 
the values of the nucleon coupling Cn' v ' T excluded by n— 208 Pb 
scattering at the 95% C. L. 



Aa M = a° xp - af = (29 ± 9) x 10~ 10 UJ, the new inter- 
action may explain this difference. From Ref. [22], 



Aa„ 



2x 2 - f3x 3 



a x 2 + (m 2 /m 2 )(l - x) 



dx . 



(13) 



where (3 = 1 for a scalar and (3 = 2 for a vector. For a ten- 
sor interaction, we trivially modify the result of Ref. [23| . 



In the limit m x <C m^, 



Cl = 47r(Aa AI /3) 1/2 < 4.8 x 10~ 4 , 
= 4 7 r(Aa AI /2) 1 / 2 < 5.9 x 10~ 4 , 
Cj = 47r(3A aAl /20) 



1/2 < 



3.2 x 10" 



(14) 
(15) 
(16) 



where the one-sided upper bounds are at the 95% C. L. 
From Fig. [5] it is evident that Eqs. (fT4")) - ([l"6"]) apply for 



m x < 10 MeV. 



The bound in Eq. (JTT]) can be combined with those in 
Eqs. (jl"4")) - (jl"6"f to give the following 95% C. L. constraints 



for m x < 10 MeV: 



a x /a < (m x /2847)^ 
a x /a < (m x /2573) 2 
a x /a < (m x /3477) 2 



scalar , 
vector , 
tensor , 



(17) 
(18) 
(19) 



with m x in MeV. A similar (numerical) procedure can be 
applied for the entire range of m x to obtain the upper 
bounds shown in Fig.[TJ We see that a vector \ with mass 
between 25 MeV (with a x ~ 10- 4 a) and 210 MeV (with 
a x ~ 10 _2 a) is a viable candidate. While a scalar \ with 
mass between 70 MeV and 210 MeV (with a x - (10~ 3 - 
10 _2 )a) is marginally allowed, a spin-2 \ is excluded. 

Muonic atom transitions. Measurements of the 
muonic 3D 5 / 2 — 2P3/2 X-ray transition in 24 Mg and 28 Si 
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atoms directly constrain a x for scalar, vector and tensor 
particles 24J. For the Yukawa form of Eq. ([5]) with the 
coupling enhanced by a factor of A, the shift in the dif- 
ference in energy levels from the QED expectation is 24| , 



AS 



T = ^ [9/(2) -4/(3)], 



(20) 



where f(j) = [1 + jm x / {2aZm^)]- 2: > , Z is the 
atomic number, and j is the principle quantum num- 
ber of the muonic state. The measured value ob- 
tained by averaging the results for 24 Mg and 28 Si, 
A£/£ = (0.2 ±3.1) ■ 10" 6 0, gives the 95% C. L. 
bound (dotted curve) in Fig. [T] No additional area of the 
relevant parameter space is excluded by this constraint. 

J/i/> decay. For m x < 2m M , the decay of scalar 
X — > is kinematically forbidden so that a constraint 

from the nonobservance of the decay J/tp — > jx-> with \ 
invisible [25} , may be employed to exclude the marginally 
allowed region with 70 MeV < m x < 210 MeV; prelim- 
inary data also exist for the decay T(35') — > 7X [26j | . 
A trivial modification of Eq. ([7|) applies to J/tp de- 
cay. The 90% C. L. upper limit on BF(J/i/j — > jx) is 
~ 4.5 x 10 -6 [2^], which when combined with BF(J/ip — > 
M + M~) = (5.93±0.06)-10- 2 [27], gives < 0.029. Then, 
the muonic Lamb shift dictates that be larger than 
3.4 x 10~ 3 which is excluded at the 95% C. L. by a p ; see 
Fig. [21 Thus, scalars are also disfavored. 

Pion decay. The 90% C. L. experimental upper limit 
on the decay ir° 7%, where x 1S a vector particle, is 
(3.3 - 1.9) • 10~ 5 for m x ranging from to 120 MeV (2^ . 
Equivalently, C% < 4.5 x 10- 4 (l-m 2 /m 2 )- 3 / 2 [U, and 
the corresponding values of CY required to explain the 
muonic Lamb shift are excluded by a M . This leaves m x 
between 120 MeV and 230 MeV. 

Eta decay. For vector x, the 90% C. L. experimental 
upper limit on invisible decays, BF(r/ — > XX)I BF(r] — > 
77) < 1.65 x 10~ 3 [30], translates into CY < 0.05 
The corresponding CY for 120 MeV< m x < m v /2 ~ 
274 MeV is excluded by a M , so that vector x is ruled out. 

Conclusions. We have found that new spin-0, spin-1 
or spin-2 particles that mediate flavor-conserving nonuni- 
versal spin-independent interactions are excluded by sev- 
eral low energy constraints as an explanation of the pro- 
ton radius anomaly. We assumed that among leptons, 
the new particles couple only to the muon so as to avoid 
the large number of constraints involving the interaction 
of the electron with exotica. We also supposed that the 
coupling of the new particle to nucleons represents the 
minimal hadronic coupling, and employed it to mesons. 

There are ways to relax some of the bounds at the ex- 
pense of introducing complication. For example, since 
the contributions of scalars and pseudoscalars to a M are 
opposite in sign, allowing both a scalar boson and a 
pseudoscalar boson with appropriately tuned couplings 
can lead to a cancellation that permits a rather large 



muonic coupling. Then, although the hadronic couplings 
are highly restricted, the muonic Lamb shift can be ac- 
commodated. Another possibility is that the new inter- 
action violates isospin or CP, so that additional freedom 
is garnered. 
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